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COMPUTER  PROGRAM 


FOR  RESERVOIR-WATER  BUDGETS 

By  Russell  R.  Schoof1 

ABSTRACT 

A  computer  program  for  computing  reservoir-water  budgets  is  pre- 
sented. The  program  was  designed  for  floodwater-retarding  reservoirs, 
but  the  system  is  applicable  to  all  water-storage  reservoirs  where  there  is 
a  need  for  computation  and  tabulation  of  water  input,  output,  and  change 
in  storage.  Discussions  of  the  computing  procedure,  input  and  output 
formats,  and  data-editing  routine  are  included.  Flow  charts  of  the  main 
program  and  two  of  the  subroutines,  complete  FORTRAN  listings,  and 
symbol  definitions  are  also  included.  KEYWORDS:  base  flow,  electronic 
computer  program,  evapotranspiration,  floodwater-retarding  reservoir, 
hydrology,  reservoir-water  budget,  runoff,  seepage,  stormflow. 


INTRODUCTION 

The  water  budget  of  a  reservoir  includes 
increments  of  change  in  reservoir  storage  and 
all  water  that  either  adds  to  or  subtracts  from 
the  storage.  Input  to  the  reservoir  includes  rain 
on  the  pool  surface  and  inflow.  Output  includes 
spillway  outflow,  evapotranspiration,  and  seep- 
age. 

Several  thousand  floodwater-retarding  struc- 
tures have  been  constructed  in  the  United 
States,  and  several  thousand  more  will  likely 
be  constructed  in  the  future.  The  Southern 
Great  Plains  Research  Watershed  was  estab- 
lished to  provide  answers  concerning  the  down- 
stream effects  of  an  upstream  flood-control 
program,  including  floodwater-retarding  reser- 
voirs. Therefore,  a  necessary  part  of  the  re- 
search was  the  development  of  a  computer 
program  for  computing  reservoir-water  budg- 
ets. 

The  program  is  written  in  FORTRAN  IV 
for  use  on  an  IBM  1130  electronic  computer 
with  8,000  words  of  core  storage  and  is  adapt- 
able to  other  computers  with  as  much  or  greater 


1  Hydraulic  engineer,  Agricultural  Research  Service, 
U.S.  Department  of  Agriculture,  Chickasha,  Okla. 
73018. 


storage.  The  program  will  be  useful  to  reservoir 
managers  and  those  engaged  in  reservoir  re- 
search. The  usefulness  of  the  program  is  not 
limited  by  the  size  of  the  reservoir. 

The  U.S.  Army  Corps  of  Engineers  at  Tulsa, 
Okla.,  computes  daily  inflow  to  its  reservoirs 
from  outflow,  change  in  storage,  and  pan  evapo- 
ration. However,  these  computations  are  done 
without  the  aid  of  an  electronic  computer.  The 
Soil  Conservation  Service  at  Stillwater,  Okla., 
and  the  Bureau  of  Reclamation  at  Amarillo, 
Tex.,  have  computer  programs  for  estimating 
monthly  water  budgets  of  reservoirs.  However, 
these  programs  are  used  only  as  a  planning 
tool  for  proposed  reservoir  sites. 

Woolhiser  and  Saxton2  and  Shanholtz  and 
Burford3  have  options  in  their  runoff-compu- 
tation programs  for  ponding  correction.  This 
option  is  applied  to  compute  instantaneous  run- 
off from  small  watersheds.  Woolhiser  and 
Saxton  also  adjust  for  rainfall  on  the  pond 


2  Woolhiser,  D.  A.,  and  Saxton,  K.  E.  1965.  Computer 
program  for  the  reduction  and  preliminary  analysis  of 
runoff  data.  U.S.  Dep.  Agric,  Agric.  Res.  Serv.  [Rep.] 
ARS  41-109,  33  pp. 

3  Shanholtz,  V.  O.,  and  Burford,  J.  B.  1967.  Computer 
systems  for  reduction  and  analysis  of  hydrologic  data. 
U.S.  Dep.  Agric,  Agric.  Res.  Serv.  [Rep.]  ARS  41-132, 
90  pp. 
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and  a  variable-drainage  area  caused  by  changes 
in  the  surface  area  of  the  pond  as  the  storage 
changes.  However,  no  other  hydrologic  com- 
puter program  has  been  found  that  computes 
daily  evapotranspiration  and  seepage  loss  from 
a  reservoir  as  a  residual  in  the  hydrologic 
equation  during  periods  of  base  inflow  or  no 
inflow.  The  program  described  in  this  report 
seems  to  be  unique  in  that  respect. 

WATER  BUDGET 
COMPUTATIONS 

If  there  is  no  base  flow  into  a  reservoir  and 
the  perimeter  of  the  pool  is  only  insignificantly 
permeable,  the  records  necessary  for  compu- 
tation of  the  water  budget  are  precipitation  on 
the  pool,  water  level  of  the  pool  (hereafter 
referred  to  as  gage  height) ,  a  gage-height  dis- 
charge rating  for  the  spillway,  and  a  gage 
height-volume  relation.  These  data  permit  com- 
putation of  the  combined  evaporation  and 
seepage  loss  from  the  pool.  If  there  is  base 
flow  into  the  pool,  it  must  also  be  measured 
and  included  in  the  computation  of  seepage  and 
evapotranspiration.  A  large  error  can  result 
in  the  water  budget  during  inflow  computation 
if  the  periphery  of  the  pool  is  extremely  porous. 
The  seepage  out  cannot  be  computed,  but  the 
return  flow  following  the  runoff  event  is  com- 
puted as  a  negative  loss.  Therefore,  if  the 
material  around  the  pool  is  highly  porous,  a 
complete  record  of  storm  inflow  as  well  as  base 
flow  is  necessary  to  determine  an  accurate 
water  budget. 

The  sediment  pools  of  floodwater-retarding 
reservoirs  are  usually  designed  to  hold  50 
years  of  expected  sediment  accumulation.  The 
pool  surface  area  and  storage  will  decrease  as 
the  reservoir  fills  with  sediment.  Therefore, 
the  storage  and  surface-area  tables  should  be 
updated  by  a  current  survey  when  a  significant 
change  in  the  tabular  values  is  suspected. 

This  program  computes  water  budgets  in 
acre-feet  with  the  hydrologic  equation  tradi- 
tionally stated  as  7=0+aS+L,  where  I  is  the 
inflow,  0  is  the  outflow,  AS  is  the  change  in 
storage,  and  L  includes  all  losses  from  storage. 
Inflow  includes  both  inflow  (Q,)  and  precipi- 
tation on  the  pool  surface  (P),  or  /=Q;+P. 

Two  procedures  are  applied  in  the  compu- 
tations. For  periods  when  the  inflow  is  base 
flow  only  or  when  there  is  no  inflow,  the  time 
increments  are  24  hours.  Then  the  water  taken 


from  storage  by  evapotranspiration  and  seepage 
(ES),  or  loss  (L)  in  the  first  equation,  is 
computed  with  Q0  substituted  for  0  in  that 
equation  to  give  the  equation  ES=I — Q0 — AS. 

During  storm  events  with  unmeasured  sur- 
face inflow,  the  program  computes  the  surface 
inflow,  but  it  must  assume  that  evaporation 
and  seepage  are  negligible  because  two  un- 
knowns cannot  be  computed  in  one  equation. 
The  data  entry  time  increments  may  range 
from  0.1  hour  to  several  hours,  depending  on 
the  rate  of  change  in  gage  height.  The  equation 
for  runoff  periods  is  Q,=Q0+aS — P.  Although 
a  negligible  ES  is  assumed  during  days  with 
inflow  events,  a  small  amount  of  ES  usually 
appears  in  the  daily  budget  printout  for  storm 
periods.  This  effect  may  be  caused  by  unequal 
time  intervals  in  the  three-point  computing 
procedure.  Also,  when  an  increment  of  inflow 
is  negative,  the  increment  is  reset  to  zero. 
Subroutine  SUMRE  then  computes  ES  for  the 
day  from  accumulated  values  of  inflow,  outflow, 
and  change  in  storage.  A  more  accurate  water 
budget  probably  results  because  the  assumption 
of  negligible  ES  during  days  with  storm  runoff 
is  incorrect. 

The  computations  during  storm  intervals  are 
made  from  the  water-level  gage  height  and 
accumulated-rainfall  entries  at  three  successive 
points.  After  completion  of  all  computations, 
associated  with  point  2,  each  point  is  advanced 
one  increment.  Thus  point  1  is  moved  to  what 
was  point  2,  point  2  to  what  was  point  3,  and- 
point  3  to  the  next  data  entry.  The  three-point 
method  was  used  to  make  the  inflow  output 
time  correspond  with  the  input  time.  The  basic 
equations  and  units  of  measurement  in  thei 
computations  are  given  in  table  1.  (The  symbols 
are  defined  in  appendix  A.) 

An  abbreviated  flow  chart  of  the  entire 
process  is  shown  in  figure  1.  Flow  charts  for 
the  main  program  and  subroutine  CHECK  are 
given  in  appendixes  B  and  C,  respectively. 

DATA  PROCESSING 

All  input  is  punched  on  a  standard  80-column 
punch  card.  The  input  includes  rating  tables 
for  reservoir  storage,  principal  spillway  out- 
flow, and  reservoir  surface  area.  Miscellaneous 
station  data  are  inserted  on  one  card.  These 
data  are  then  followed  by  a  station-year  of 
hydrologic  data  consisting  of  reservoir  gage 
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Figure  1. — Flow  chart  of  water  budget  computation  process. 
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height,  accumulated  rainfall,  base  inflow,  and 
the  associated  dates  and  times. 

A  preliminary  edit  of  each  year  of  input  data 
is  made  by  the  EDIT  program  listed  in  ap- 
pendix D.  This  program  examines  each  data 
entry.  Error  messages  are  printed  for  days 
with  data  missing,  cards  out  of  sequence,  and 
gage-height  increments  greater  than  0.4  foot 
between  entries.  The  amount  of  rainfall  re- 
ceived each  day  is  computed,  stored,  and  then 
printed  on  one  sheet  where  it  can  be  easily 
compared  with  tabulated  daily  rainfall.  Sub- 
routine CHECK  (appendixes  C  and  F)  edits 
the  data  again  just  before  the  water  budget 
computations  are  made.  The  second  editing  is 
to  ascertain  that  no  cards  have  been  lost  and 
that  all  errors  in  the  data  have  been  satis- 
factorily corrected.  Subroutine  CHECK  does 
everything  the  EDIT  program  does  except 
compute  and  tabulate  the  daily  rainfall. 

Station  ratings  and  data  errors  are  printed 
by  subroutine  CHECK.  Storm-event  output  data 
are  printed  by  the  reservoir-water  budget  main 
program.  Monthly  and  annual  summaries  are 
printed  by  subroutine  SUMRE.  Although  only 
printed  output  is  given,  the  programs  could  be 
easily  modified  to  include  output  in  any  other 
mode  desired. 

Input 

Rating  tables 

Storage,  outflow,  and  surface  area  are  com- 
puted with  a  parabolic  interpolation  subroutine, 
SPPNT  (see  appendix  G),  developed  by  Du- 


Chateau  et  al.4  The  routine  uses  six  consecutive 
values   (three  on  each  side  of  the  point  in  j 
question)  from  the  appropriate  table.  There-, 
fore,  each  table  must  begin  with  three  non- 
positive  values.  The  third  entry  is  zero.  Each  j 
rating-table  card  contains  eight  values  with 
two  places  right  of  the  decimal  (format  8F8.2). 
The  last  entry  in  the  table  must  be  blank  or 
zero  to  signify  the  end  of  the  table.  Therefore, 
if  the  number  of  tabular  entries  is  a  multiple'1 
of  eight,  a  blank  card  must  be  included  at  thej 
end.  Columns  65  through  80  can  be  used  for,) 
the  table  identification  and  card  sequence,  but II 
they  are  not  read  by  the  computer. 

All  rating-table  input  is  punched  on  the 
standard  80-column   punch  card   in   format] ! 
8F8.2.  The  following  data  are  included: 

1.  Reservoir-storage  rating  table  in  acre- 
feet  for  gage-height  increments  of  1.0  foot. 

2.  Principal  spillway  outflow  rating  in  cubic 
feet  per  second  for  0.2-foot  gage-height  incre- 1 
ments. 

3.  Reservoir-surface-area  table  in  acres  foril 
1.0-foot  gage-height  increments  corresponding 
to  gage-height  points  in  the  storage  ratingii 
table. 

Station  data 

The  first  three  values  on  the  station  data  J 
card  are  data  identifiers.  VMIN  is  subtracted, 
from  the  gage-height  input  values  to  obtain  a] 
value  that  is  then  used  to  select  the  correct 


4  DuChateau,  P.  C;  Nofziger,  D.  L. ;  Ahuja,  L.  R. 
and  Swartzendruber,  D.  1972.  Experimental  curves  and:  j 
rates  of  change  from  piecewise  parabolic  fits.  Agron.  J.  I 
64:  538-542. 


Table  1. — Measurements  and  equations  for  water  budget  computations 
Measurement  Equation 

Time  lapse  between  points  1  and  3   h  ...  DT  =  T3-T1. 

Time  lapse  between  points  1  and  2   h     .  DELT  =  T2  — Tl. 

Average  outflow  rate  from  point  1  to  point  2  ■■■■  ft3/s    ; .  QPA=  (QP1  +  QP2)  /2.O. 

Volume  of  outflow  between  points  1  and  2   acre  f t  ■  •  ■  DQP  =  QPA*DELT/12.1. 

Total  change  in  storage  from  point  1  to  point  3  .  acre  f t  ■  ■  •  DV  =  VAF3  —  VAF1. 

Rainfall  on  pool  between  points  1  and  3   acre  ft  •  ■•  PV=  (PA3-PA1)  *A2/12.0. 

Rate  of  inflow  assigned  to  point  2   ftVs  ■  •  •  QI2  =  QP2+ [(DV-PV)  / 

DT]*12.1. 

Volume  of  inflow  between  points  1  and  2   ft3/s-h  •■•  DI=  [(QI2  + QI1)  / 

2.0]*DELT. 

Accumulated  inflow  volume  to  point  2   ft3/s-h  ...  ADV  =  ADV  +  DI. 

Watershed  inflow  between  points  1  and  2   in       DI=  (DI/12.1)  *12.0/ 

(AREA  — A2). 

Accumulated  watershed  inflow  to  point  2   in  •  • .  AII  =  AII  +  DI. 
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surface-area  and  reservoir-storage  values  for 
the  parabolic  interpolation  subroutine,  SPPNT. 
QMIN  has  a  similar  function  in  the  computa- 
tion of  spillway  outflow.  BP  is  compared  with 
each  gage  height  to  determine  if  there  is 
principal  spillway  flow.  BE  is  compared  with 
each  gage-height  value  greater  than  BP  to 
determine  if  there  is  flow  through  the  emer- 
gency spillway.  This  program  does  not  include 
a  provision  for  computation  of  flow  through 
emergency  spillways.  However,  with  a  slight 
modification  in  the  program,  emergency-spill- 
way flow  could  be  included  in  the  principal 
spillway  rating. 

Watershed  area  (A)  is  used  in  converting 
inflow  in  cubic  feet  per  second  to  water  yield 
in  inches.  The  pool  surface  area  is  subtracted 
from  A  before  making  each  computation  of 
water  yield.  The  Julian  calendar  day  (NUM) 
is  the  initial  point  for  a  day  counter  (1  if 
processing  starts  on  January  1). 

The  station  data  card  is  prepared  as  shown 
in  table  2. 

Hydrologic  data 

The  first  three  values  of  each  hydrologic 
data  card  are  for  identification  only  and  are 
not  read  into  the  computer.  For  days  with  in- 
sufficient precipitation  to  produce  runoff  and 
less  than  0.4  foot  of  change  in  gage  height, 
entries  are  input  at  midnight  only.  However, 
the  time  is  designated  at  25.0  to  distinguish 
these  days  from  midnight  entries  for  days  with 
storm  inflow.  The  precipitation  is  the  accumu- 

Table  2. — Card  format  for  station  data  input 


Data  description 

Symbol 

Unit 

Format 

Reservoir  No.   

NR 

13 

Watershed  No.   

NWS 

14 

Year   

NY 

13 

Gage  height  corresponding  to 

first  value  in  storage  and 

surface  area  tables  less 

1.0  foot  

VMIN 

Feet 

F6.2 

Gage  height  corresponding  to 

first  entry  in  principal 

spillway  rating  less  0-2  foot 

QMIN 

Feet 

F6.2 

Gage  height  of  principal 

spillway  

BP 

Feet 

F6.2 

Gage  height  of  emergency 

spillway  

BE 

Feet 

F6.2 

Watershed  area  at  dam   

A 

Acres 

F8.1 

Julian  calendar  day  No.  of 

first  day  of  data  input  

NUM 

14 

lated  amount  during  the  previous  24  hours.  The 
base  flow  punched  on  these  25.0  time  cards  is 
the  average  base  flow  from  all  combined  inflow 
channels  during  the  previous  24  hours.  Base 
inflow  is  input  as  zero  or  is  left  blank  on  days 
with  storm  inflow. 

Days  with  storm  inflow  following  a  single- 
card  input  day  begin  with  the  25.0  time  card 
for  the  previous  day.  However,  subdivided  days 
following  a  subdivided  day  must  begin  with 
0.0  time  input.  The  rainfall  input  on  the  last 
single-card  input  day  preceding  each  storm 
event  and  on  subdivided  days  is  an  accumulated 
amount.  The  accumulation  can  restart  from 
zero  at  any  time.  However,  the  first  value  in  a 
new  accumulation  must  be  smaller  than  the 
final  value  of  the  preceding  accumulation. 
Usually,  the  rainfall  accumulation  is  restarted 
at  the  beginning  of  each  rainfall  chart. 

The  card  format  for  the  hydrologic  data  input 
is  given  in  table  3. 


End  of  data 


Hydrologic  data 
for  1  year 


Station  data  card 
_/  Reservoir  surface  area  table 

Reservoir  storage  rating  table 

Figure  2. — Proper  stacking  of  input  cards  for  each 
station-year. 


Table  3. — Card  format  for  hydrologic  data 
input 


Data  description 


Reservoir  No. 
Watershed  No. 

Year   

Month   

Day   


Symbol 

Unit 

Format 

NR 

13 

NWS 

14 

NY 

13 

NMO 

13 

NDAY 

13 

Tl,  T2,or  T3 

Hours 

F5.1 

Gl,  G2,  or  G3 

Feet 

F6.2 

PI,  P2,  or  P3 

Inches 

F6.2 

Accumulated 

rainfall   

Base  inflow   BF1,  BF2,  or  BF3  ftVs  F6.2 
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VOLUME  RATING  AT 

ONE  FOOT 

INTERVALS   IN  ACRE  FEET 

-1.00 

-0.50 

0.00 

0.20 

0.50 

1.00 

2.00 

4*00 

7.00 

10.50 

15.00 

20.00 

26.00 

32.00 

40  .  00 

48.UU 

58.00 

70.00 

83.00 

98  .00 

113.00 

130.00 

l  c. n  .  oo 

17  3.00 

197.00 

225.00 

256.00 

290.00 

328.00 

370.00 

n .  no 

PRINCIPAL 

SPILLWAY  RATING 

AT  0.2 

FOOT   INTERVALS  IN 

r  c  c 

-0.08 

-0.04 

0.00 

U«25 

0.93 

1  •  97 

3  .  1  © 

5.30 

6.15 

7.20 

7  •  85 

8*65 

9.40 

10.10 

10.70 

11.30 

12.00 

12.55 

13  .20 

13.75 

14.40 

1/.  .OS 

16.05 

16.60 

17.15 

17.70 

18.25 

18.75 

ig.ifl 
X  Tr  .  a  U 

20.40 

20.95 

21.50 

22  .05 

22.60 

23.15 

A  J  .  f  3 

24.80 

25.40 

25.95 

26*50 

27.05 

27.65 

5  ft  .  3  =, 
ft  O  .  £  3 

O.  OO 

SURFACE  AREA  RATING  AT  ONE 

FOOT 

INTERVALS 

IN  ACRES 

-0.30 

-0.10 

0.00 

0.10 

0.30 

0.70 

1  .  5  U 

3  .  ^.0 

3.30 

4*00 

4.80 

5*30 

6.00 

6.80 

8.00 

9.90 

10.90 

12.10 

13.20 

14.60 

16.60 

19.10 

21.10 

23.00 

25.80 

29.60 

32.90 

35.60 

38.70 

42.00 

u  •  uu 

GAGE 

HEIGHT 

INCREMENT 

GREATER  THAN  0. 

40  FEET 

7  621  75 

5  22 

19.7  20. 

43     2.64  0.00 
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Figure  3. — Output  from  subroutine  CHECK. 
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Figure  4. — Output  for  each  period  of  1  or  more  days  of  surface  inflow. 


If  December  31  has  a  single-card  input  with 
:25.0  time,  it  is  followed  by  a  blank  card.  How- 
ever, if  December  31  is  a  subdivided  day,  the 
tnext  to  last  card  is  an  entry  for  0.0  time  on 
January  1,  and  the  last  card  is  an  entry  for 
the  next  data  point  with  a  zero  or  blank  in 
the  month  space  to  signify  the  end  of  data. 
The  input  cards  for  each  station-year  of  data 
must  be  stacked  as  shown  in  figure  2. 

Output 

The  output  for  subroutine  CHECK  (figure 
3)  includes  listings  of  the  reservoir  volume, 
discharge  and  surface-area  ratings,  gage-height 
and  time-sequence  errors,  and  Julian  day  num- 
ber for  December  31.  If  any  errors  are  de- 
tected, the  processing  will  terminate  upon 
completion  of  the  data  checking.  The  first  line 
of  output  from  the  main  program  includes 


VMIN,  QMIN,  BP,  BE,  AREA,  and  NUM  (see 
appendix  A). 

Output  for  each  period  of  1  or  more  days  of 
surface  inflow  is  then  listed  as  shown  in  figure 
4.  The  listing  includes  the  reservoir  arid  water- 
shed numbers,  year,  month,  day,  time,  gage 
height,  rainfall  on  the  pool,  volume  of  water 
stored  in  the  pool,  rate  of  flow  from  the 
principal  spillway,  accumulated  principal-spill- 
way discharge,  rate  of  inflow,  and  accumulated 
inflow.  Although  columns  are  provided  for 
emergency-spillway  flow,  the  program  does  not 
compute  it. 

The  daily  summaries  are  stored  on  disk.  At 
the  end  of  a  station-year  of  computation,  sub- 
routine SUMRE  (see  appendix  H)  is  called  to 
complete  the  computation  of  daily,  monthly, 
and  annual  summaries  and  list  the  output  as 
shown  in  figures  5  and  6. 
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APPENDIX  A. 


SYMBOL  DEFINITIONS  FOR  WATER  BUDGET 
COMPUTATIONS 


A 

GQ1,  GQ2,  and  GQ3  divided  by  0.2. 

NDAY, 

Al, 

MDAY,  IDAY 

Day. 

A2,  A3 

Surface  area  of  reservoir  at  Tl,  T2, 

NER 

Line  counter. 

and  T3,  respectively. 

NF 

Month. 

ACP 

Accumulated  daily  precipitation 

NMO 

(inches) . 

MMO,  IMO 

Month. 

ACPS 

Accumulated  principal  spillway 

NR 

Reservoir  number. 

outflow  (acre-feet). 

NUM 

Julian  calendar  day  counter. 

ADV 

Accumulated  daily  inflow 

NWS 

Watershed  number. 

(cubic  feet  per  second-hours  and 

NY 

Year. 

acre-feet) . 

PA1, 

All 

Accumulated  inflow  (inches). 

PA2,  PA3 

Accumulated  precipitation  (inches) 

APV 

Precipitation  on  reservoir 

PI, 

(acre-feet) . 

P2,  P3 

Accumulated  precipitation  (inches) 

AQPT 

Accumulated  daily  principal-spillway 

PI 

Precipitation  increment  (inches) 

outflow  (acre-feet). 

between  Tl  and  T2. 

AREA 

Drainage  area  at  dam  including  area 

PV 

Volume  of  precipitation  falling  on 

of  pool  (acres) . 

reservoir  (acre-feet). 

BE 

Gage  height  of  emergency  spillway 

QI 

Rate  of  inflow  to  reservoir 

(feet) . 

(cubic  feet  per  second) . 

BF1, 

QI1, 

BF2,  BF3 


BP 

DELT 
DI 

DQP 

DT 

DV 


G 

GHS 

Gl, 

G2,  G3 
GQ1, 

GQ2,  GQ3 
GV1, 

GV2,  GV3 
JAB 
JD 

JOB 

M 

ND 


Base  inflow  (cubic  feet  per  second) 

at  Tl,  T2,  and  T3,  respectively. 

These  values  are  not  entered  on 

subdivided  days. 
Gage  height  of  principal  spillway 

(feet). 
T2-T1. 

Inflow  from  Tl  to  T2  (cubic  feet 
per  second-hours  and  inches) . 

Principal-spillway  outflow  from  Tl 
to  T2  (acre-feet). 

T3-T1. 

Incremental  change  in  reservoir 
contents  between  Tl  and  T3 
(acre-feet) . 

G2  + 0.001. 

Gage  height  at  start  of  computation 
period. 

Gage  height  (feet). 

Gl,  G2,  or  G3-QMIN. 

Gl,  G2,  or  G3-VMIN. 
Input  data  (file  1)  counter. 
Internal  variable  controlling 

program  branching. 
Daily  summary  data  (file  2)  counter. 
Day. 

Point  in  reservoir  volume,  surface 
area,  or  spillway  discharge  array 
that  is  at  the  left  of  the  interval  in 
which  interpolation  will  be  made. 


QI2,  QI3 
QMIN 
QPA 

QP1,  QP2 

Tl, 
T2,  T3 

VAF1, 

VAF2,  VAF3 

VMIN 

XAD 

YA 

YQ 

YV 
Z 

ZA 


Rate  of  inflow  (cubic  feet  per 

second)  at  Tl,  T2,  and  T3, 

respectively. 
0.2  foot  less  than  gage  height  of  zero 

entry  in  principal-spillway 

discharge  table. 
Average  principal-spillway  outflow 

from  Tl  to  T2. 
Rate  of  principal-spillway  outflow 

(cubic  feet  per  second). 

Time  (hours)  at  data  entries 
1,  2,  and  3. 

Reservoir  storage  (acre-feet)  at 
Tl,  T2,  and  T3,  respectively. 

1  foot  less  than  gage  height  of  zero 
entry  in  storage  and  area  rating 
tables. 

The  fractional  distance  to  where 

interpolation  is  desired. 
Reservoir  surface  area  array  (area), 

1-foot  intervals. 
Principal  spillway  discharge  array 

(cubic  feet  per  second),  0.2-foot 

intervals. 
Reservoir  volume  (acre-feet),  1-foot 

intervals. 
Array  of  mean  daily  reservoir  inflow 

(cubic  feet  per  second) . 
Continuation  or  termination  code  in 

subroutine  CHECK. 
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APPENDIX  B.-FLOW  CHART  FOR  RESERVOIR  WATER  BUDGET 


CALL  CHECK  TO 
STORE  AND  CHECK  INPUT 


WRITE  RATINGS 
AND  COLUMN  HEADINGS 


COMPUTE  POINT  1  STORAGE 

PA1  -  PI 

GHS  =  Gl 

GV1  =  Gl-VMIN 

NO  =  GV1 

XAD  =  GVl-FLOAT(IFIX(GVI)) 

CALL  SPPNT(YV,ND,XAD,VAF1 

READ 
NDAY 


(l'JAB)  NMO,  \^ 
,  T1,G1,P1,BF1 


ZERO  ACCUMULATORS 

AND 

MEAN  DAILY  INFLOW 

ARRAY 

JAB  = 

1 

JOB  = 

1 

NER  = 

0 

All  = 

.0 

ACP  = 

.0 

ACPS  = 

.0 

APV  = 

.0 

ADV  = 

.0 

AQPT  -• 

.0 

SET  Z  FII 

t:  - 

.0 

COMPUTE  POINT  1  OUTFLOW 


=  Gl-QMIN 
=  GQ1/.2 
=  A 

=  A-FLOAT(IFIX(A)) 


CALL  SPPNT(YQ,ND,XAD,QP1) 


READ  (l'JAB)  NMO 
NDAY,T2,G2,P2,BF2 


MMO 

=  NMO 

MDAY 

=  NDAY 

GV2 

=  G2-VMIN 

COMPUTE  POINT  2  STORAGE 
AND  SURFACE  AREA 

ND  =  GV2 
XAD  =  GV2-FLOAT(IFIX(GV2)) 
CALL  SPPNT(YV,ND,XAD,VAF2) 
CALL  SPPNT(YA,ND,XAD,A2) 


COMPUTE  POINT 

GQ2  =  G2-QMIN 
A  =  GQ2/.2 
ND  =  A 
XAD  =  A-FLOAT  (IFTX(A)' 
CALL  SPPNT(YQ,ND,XAD,QP2) 


AQPT 
ACPS 
ADV 
APV 
ACP 


SUMMARY 

((QP1+QP2)/12. 1)*12 
ACPS+AQPT 
BF2*2A. 
P2*A2/12. 
P2 


ADVANCE  POINT 

1  TO  POINT  2 

Gl  = 

G2 

BF1  = 

BF2 

QP1  = 

QP2 

VAF1  = 

VAF2 

PA1  = 

P2 

JD  - 

1 
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CHANGE  T2  TO  24 

0 

AND  MOVE  DAY  BACK 

ONE 

T2  =  24.0 

NMO  =  MMO 

NDAY  =  MDAY 

1MO  =  IMO 

MDAY  =  IDAY 

COMPUTE  POINT  3  STORAGE 
AND  SURFACE  AREA 

GV3  =  G3-VMIN 

ND  =  GV3 
XAD  =  GV3-FLOAT(IFIX(GV3)) 
CALL  SPPNT(YV,ND,XAD,VAF3) 
CALL  SPPNT(YA,ND,XAD,A3) 


DV  -  VAF3-VAF1 

?V  -  (PA3-PA1)*A2/12.0 

QI  -  ((DV-PV)/DT)*12.1 


WRITE  PROGRAM  DOES  NOT 
COMPUTE  EMERGENCY  SPILLWAY 
FLOW 


WRITE  POINT  2  OUTPUT 
WRITE  (3,92)  NR , NWS , NY , NMO , NDAY 
T2,G2,PI,VAF2,QP2,ACPS,qi2,AII 


DI 

((QI2+QI1)/2.)*DELT 

ADV 

ADV+DI 

DI 

(DI/12.1)*12.0/(AREA-A2) 

All 

AII+DI 

QI2 

QI2+.005 

G2 

G2+.001 

PI 

PI+.001 

QP2 

QP2+.005 

ADVANCE  POINT  1  TO  POINT  2 

AND  ACCUMULATE  STORM 

RAINFALL 

NER  -  NER+1 

Gl  -  G2-.001 

BF1  =•  BF2 

QP1  «  QP2-.005 

VAF1  -  VAF2 

Tl  -  T2 

PA1  -  PA2 

QI1  -  QI2-.005 

JD  =  2 

>  1 

200 
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APPENDIX  C.-FLOW  CHART  FOR  SUBROUTINE  CHECK 


Checks  Input  for  days  with  surface  inflow 
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APPENDIX  I).— FORTRAN  LISTING  FOR  EDIT  PROGRAM 


DIMENSION  Z (  12*31 >  »T  (  12  ) 

1  FORMAT ( I3»I4»3I3»F5.1*2F6«2*F6.2) 

2  FORMAT (1H   I  3 » 1 4 • 3  I  3 » F5 . 1 » 2 F6 . 2 »F6 . 2 ) 

3  FORMAT (  1H05X. 'GAGE   HEIGHT    INCREMENT   GREATER   THAN  0.25  FEET') 

4  FORMAT ( 1H05X » 45HCARDS  MISSING  OR   DATE-TIME  GROUP  OUT  OF  ORDER) 

5  FORMAT ( 1H1 ) 

6  FORMAT ( 1H05X* 'GH  RISE    EXCEEDS   PRECIPITATION*    INSTANTANEOUS  VALUES 
1NEEDED 1 ) 

7  FORMAT ( 1H0 ) 

8  FORMAT ( 1H   I4»1X,'DAYS    IN   THIS  DcCK') 
ICO  =  0 

PA=.  0 

DO  20  NF=1*12 
DO  2  0  M=l»31 
2C  Z(NF»M)=.0 
WRITE  (3*5) 
GO  TO  22 

21  WRITE(3»2)    LS *LWS»LYR »LMOf LDAY * T I M » G » P ♦  QI 
WRITE    (3»2)  NS»NWS*NYRtNM0»NDAY»TIME*G2»PE*QI 
WRITE (3*7  J 

22  READ    (2*1)    LS  ♦  LWS » LYR » LMO * LDAY * T I M * G 1 ♦ P » Q I 

25  I  F ( T  IM-24.C > 55*42  ,26 

26  Z ( LMO ♦ LDAY ) =P 
GO   TO   2  9 

2  7  LMO=NMO 

LDAY= 1 
GO  TO  29 

28  LDAY=NDAY 
T I M* TIME 

29  ICO=ICO+l 
P1  =  P 

READ ( 2  » 1 )NS  *NWS*NYR*NMO»NDAY*TIME »G2  *P  *QI 
PE  =  P 

IF(NS-LS)43*32*43 
32  DIF=G2-G1 
G  =  G1 
G1  =  G2 

IF(ABS(D IF) -.25) 34*34*33 

3  3  WRITE (3*3) 

WRITE (3*2)    LS »LWS*LYR »LMO»LDAY *T IM»G»P ♦  QI 

WRITE    (3*2)  NS»NWS»NYR»NMO»NDAY»TIME»G2»PE»QI 

WRITE (3*7) 
34   IF(TIME-24. 0)67,42*35 
3  5   Z ( NMO ♦ NDAY ) =P 

IF(P)"38*38  *36 

36  DIF=(DIF-.01)*12« 
IF(DIF-P)38*38*37 

37  WPITE( 3*6) 
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WRITE ( 3 t2 )  LS  »LWS»LYR  »LMO»LDAY  »TlM»GtPl »QI 
WRITE(3»2)NS»NWStNYR»NMO»NDAY»TIMEtG2»P»OI 

38  IF(NDAY-LDAY-1)40»39.42 

39  IF(NM0-LM0)42 »28»42 

40  IF(NMO-LMO-l >42f4l »42 

41  IF    (NDAY-1)  42f27»42 

42  WRITE  (3»4) 
GO  TO  21 

43  WR I TF ( 3  »  8 ) I  CO 
ICO  =  0 

WRITE  (3»5) 
9  FORMAT    (  1H  .»12F10.2  ) 

10  FORMAT (1H  F10.2»10X»10F10.2) 

11  FORMAT    (1H   »F10.2 »10X »F10.2 » 10X » F 10 • 2 . 1 OX ♦ F 1 0 . 2 »F10.2 »10X »F10«2 »10 
1X»F10.2) 

12  FORMAT ( 1H030X , 'DAILY   PRECIPITATION   FOR  SITE'»I4»«    WAT ERSHED 1  , 1 4 . • 
1YFAR ' t 14 ) 

13  FORMAT    (1H0120H  JAN  FEB  MAR  APR  MAY 
1           JUNE             JULY                AUG             SEPT                OCT  NOV 
2DEC) 

WRITE(3»12)LS»LWS»LYR 

WRITE(3»13) 

M  =  l 

NF=1 

44  <  =  0 

45  WRITE    (3»9) (Z(NF»M) #NF=1»12) 
M  =  M+  1 

K  =  K+1 

IF    (  N'-29)46»48»46 

46  IF    ( K-5 ) 45 »47 f45 
4?  W  R I T  E ( 3  *  7 ) 

GO   TO  44 

48  'a'  R  I  T  E    (3»10)Z(1»29)  ♦  (Z(NF»M)  »NF=3tl2) 
M  =  M+1 

WRITE    (3»10)Z  (1.30)  »_(Z(NF»M)  tNF=3»12) 
M  =  M+1 

WRITE  (3»ll)Z(l,31)»Z(3»31)»Z(5.31)»Z(7t31),Z(8»31)tZ(10»31)»Z(12 
1  »31  ) 
SUM1=.C 
00  50  NF=1,12 
DO  49  M=1.31 

49  SUM  1  =  SUM1  +  Z ( NF  »M ) 
T  ( ,MF  )  =SUM1 

50  SUM1=.0 

WRITE    (3*9) (T(NF) »NF=1,12) 

WRITE( 3»5) 

LWS=NWS 

LYR=NYR 

LS  =  NS 

LMC=NMO 

LDAY=NDAY 

T  I  M=  T I  ME 

G1=G2 

DO  51  NF=lil2 
DO  51  M=lf31 

51  Z(NF»M)=.0 
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GO  TO  23 

52  LMO=NMO 
LDAY= 1 
ICO=lCO  +1 
GO  TO  54 

53  LDAY=NDAY 
ICO= ICO+1 

54  T I M  =  T I  ME 
P1  =  P2 

P  =  P2 

55  READ    (2*1)     NS • NWS » NYR ♦ NMO t NDAY » T I ME *G2 » P2 »Q I 

PE  =  P2 

IF(NWS-LWS)43 »56»43 

56  CALL   DATSW  (1»J) 
IF1J-1!  57*58*57 

57  WRITE    I  3  » 2  >  NS  *NWS » NYR ♦ NMO » NDAY  » T I  ME  » G2  » P  2 ♦ Q I 

58  DIF=G1-G2 
G  =  G1 

G1  =  G2 

IF(ABS(D IF)-. 25)60*60*59 

59  WRITE(3»3) 

WRITE(3»2)    LS»LWS»LYR»LM0»LDAY»TIM»G»P»  QI 
WRITE    (3»2)   NS  *NWS  »NYR  t NMO  t NDAY ♦TIME»G2»PE»QI 
WRITE (3*7) 

60  IF( TIME-24.0) 61*42 »75 

61  IF(P2-P1)63. 62*62 

62  PA=PA+P2-P1 
GO  TO  64 

63  PA=PA+P2 

64  IF(TIME-TIM)71*42»65 

65  IF(NDAY-LDAY)42»66»42 

66  I  F(NMO-LMO) 42*54*42 

67  ICO=ICO-l 
IF(P-P1 )68»69»69 

68  PA=P 

GO  TO  70 

69  PA=P-P1 

70  P2=P 
MMO=  NMO 
MDAY=NDAY 
GO  TO  73 

71  IF(TIME)42»72*42 

72  Z ( MMC  »  MDAY ) =P A 
PA=*0 

MMO=  NMO 
MDAY=NDAY 

73  IF(NDAY-LDAY-1)74*78.42 

74  IF(NMO-LMO-l )^2*52 »42 

75  IF(NDAY-LDAY)42*76*42 

76  IF(TIM)42*77*42 

77  Z ( NMO  »  NDAY ) =P 2 
P  =  P2 

GO  TO  39 

78  IF(NMO-LMO)42*53*42 

END 
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APPENDIX  E.-FORTRAN  LISTING  FOR  RESERVOIR- WATER 

BUDGET 


DEFINE   FILEU  1280  » 10  »U  .  JAB  ) 
DEFINE   FILE2( 549 ♦ 1 5 » U » JOB > 

DIMENSION   T(12)  tZ(  12*31)  tYVUO)  »YQ(  100)  »YA(40) 
COMMON   JAB  *  JOB 
79  FORMA T ( 1H  ) 

83  FORMAT ( 1H  4F7  .2  ,  F7 . 1 , I  4 ) 

84  FORMAT ( lH    1    NR  STA   YR   MO   DY   TIME        GH     RAIN        VOL  PR   SP  ACPSPD 
1   EM   SP     ACESPD      IN  FL  ACCINFL') 

85  FORMA T ( lH    '  (HR)      (FT)    (IN)      (ACFT)       (CFS)  (ACFT) 
1    (CFS)      (ACFT)    (CFS)  (IN)') 

90   FORMAT (1H    'ROUTINE   FOR   EMERGENCY   SPILLWAY   FLOW  NOT    IN   PROGRAM  1  ) 
92   FORMAT ( lH   1 3 f 1 4 » 3  I  3 » F5 • 1 » F 6« 2 • F5 • 2 • F8 . 2  » F6 • 2  * F8 .2  » 14X  »F7 • 2 » F8 • 5 ) 

C  READ  AND   STORE   DATA  ON   DISK   AND  CHECK  DATA 

95   CALL   CHECK    ( Y V ♦ YQ » NR ♦ N WS * N Y » VM I N » QM I N » BP ♦ BE » ARE A » NUM » YA » GHS ) 


C  WRITE   HEADER   CARD   DATA   AND   COLUMN  HEADINGS 
WRITE! 3  1 83 ) VM IN »QM I N *   BP • BE » AREA t NUM 
WRITE! 3»S4) 
WRITE (3*35) 

C  ZERO  ACCUMULATORS  AND  MEAN   DAILY    INFLOW  ARRAY 
JAB  =  1 
JOB=l 
N  E  R  =  0 
A  I  I  =  .  0 
ACP=.0 
ACPS=.0 
APV=.0 
ADV= .0 
AQPT=.C 

DO   100  NF=1»12 
DO   100  M=l,31 
IOC  Z(NF.M)=.0 

C  READ   POINT    1  INPUT 

R EAD ( 1 1  JAB ) NMO*NDAY*  Tl »G1 ♦  PI tBFl 

C  COMPUTE   POINT   1   STORAGE  AND  OUTFLOW 
PA1=P1 
GHS=G1 
GV1=G1-VMIN 
ND=GV1 

X AD  =  GV 1-FLOAT ( IF  IX ( GV1  )  ) 
CALL   SPPNTt YV »ND.XAD»VAF1  ) 
IF(Gl-BP) 109»109.108 
108  GC1=G1-QMIN 
A=G01/.2 


21 


N  D  =  A 

XAD= A-FLOAT ( I F I X ( A )  ) 
CALL  SPPNT ( YQ  »ND  »X AD  »QP1  ) 
GO  TO  110 

109  OP1=.0 

C   PEAD  POINT   2  INPUT 

110  READ ( 1 1  JAB) NM0*NDAY»T2 *G2  *P2  *BF2 
MMO  =  NMO 

MDAY=NDAY 

GV2=G2-VM1N 

I F ( NMO ) 114  »400t 1 14 

C  COMPUTE   POINT   2   STORAGE  AND  SURFACE  AREA 
114  ND=GV2 

XAD  =  GV2-FLOAT ( IF  IX (GV2  )  ) 
CALL   SPPNT ( YV »ND»XAD tVAF2 ) 
CALL   SPPNT    { Y A  *  ND  *  XAD  *  A2 ) 
IF(T2-24.)12l»l2l»116 

C   COMPUTE   POINT  2  OUTFLOW 

116  IF(G2-BP)U7»117.118 

117  QP2=.0 

GO  TO  120 

118  GQ2=G2-QMIN 
A  =  GQ2 /  •  2 

ND  =  A 

XAD  =  A-FLOAT ( I  FIX ( A )  ) 
CALL  SPPNT ( YQ »ND*XAD »OP2 ) 

C   COMPUTE   DAILY   SUMMARY   FOR   DAYS   WITH  MIDNIGHT   POINTS  ONLY 

120  AOPT= ( (QP1+GP2 ) /12. 1 )*12. 
ACPS=ACPS+AOPT 
ADV=BF2*24. 
APV=P2*A2/12. 

ACP=P2 

C  ADVANCE   POINT   1   TO  POINT  2 
G1  =  G2 
BF1=BF2 
QP1=QP2 
VAF1=VAF2 
PA1=P2 
JD  =  1 

GO  TO  200 

C   COMPUTATION   FOR   DAYS  WITH   SURFACE  INFLOW 

C   ZERO   DAILY  ACCUMULATORS 

121  £CP=.0 
APV= .0 
AOPT=.0 
T1  =  .0 
QI  1=.0 
PA2=P2 
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C  WRITE   COLUMN  HEADINGS 
WRITE! 3*79) 
IF(NER-50>123.123.122 

122  WRITE13.84) 
WRITE (3*85) 
NER  =  0 

123  IFIP2-PA1 )124.l33.133 

124  PA1=.0 

C  READ  POINT   3  INPUT 

133  R  EAD ( 1 • JAB )    I  MO , I  DA Y ♦ T 3 . G3 , P3 . BF 3 

C  COMPUTE   POINT   3   STORM   ACCUMULATED  PRECIPITATION 
IFIP3-P2) 137.133.138 

137  PA3=PA2+P3 
GO  TO  141 

138  PA3=PA2+P3-P2 

HI    IF< T2  )  143.142  .143 

C   IF   T2    IS    .0   CHANGE   TO   24.0   AND   MOVE   DAY   BACK  ONE 

142  T2=24. 

NMO=MMO 
NDA Y=MDAY 
MMO= I  MO 
MDAY= I  DAY 

C  COMPUTE   POINT   3   STORAGE   AND   SURFACE  AREA 

143  GV3=G3-VMIN 
ND=GV3 

XAD  =  CV3-FLOAT ( IF  IX (GV3  )  ) 
CALL   SPPNT    ( YV.ND.XAD.VAF3  ) 
CALL   SPPNT    ( YA.ND.XAD.A3 ) 

C  COMMUTE    INST.    RATE  OF   CHANGE    IN   STORAGE   AT   PT.    2    EXCLUDING  PRECIP.    ON  POOL 


I  F 

(T3)  151.150.151 

150 

DT  = 

24. -Tl 

GO 

TO  157 

151 

IF 

(T2-24.0)  155.154.154 

154 

DT  = 

DEL  T  +  T3 

GO 

TC  157 

155 

DT  = 

T3-T1 

157 

DV= 

VAF  3-VAF 1 

PV  = 

( PA3-FA1 )*A2/12. 

01  = 

(  ( DV-PV ) /DT  )  *  1 2  •  1 

C  COMPUTE   AND  ACCUM.   PRECIP.    ON  POOL 
PI=(PA2-PA1 ) 
ACP=ACP+PI 
PV=PI*A2/12. 
APV=APV+PV 

IF    (G2-BP)  162.162.163 

162  QP2=.0 

GO  TC  167 

163  IF    (02-BE)  166.166.164 
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C    IF   THERE    IS   EMERGENCY   SPILLWAY   FLOW »   PRINT   MESSAGE  AND  BYPASS 
C   REMAINDER   OF   DATA  COMPUTATION 
164  WRITE  (3»90) 
GO   TO  400 

C   COMPUTE   AND  ACCUMULATE   POINT   2  OUTFLOW 

166  GQ2  =  G2-QM I N 
A=GQ2/.2 

N  D  =  A 

XAD  =  A-FLOAT ( I  FIX (A )  ) 
CALL  SPPNT ( YQ »ND»XAD»QP2 ) 

167  QPA= (OP1+QP2 ) /2.C 
DQP=QPA*DELT/ 12. 1 
AGPT=AQPT+DQP 
ACPS=ACPS+DGP 

C   COMPUTE   POINT   2    INFLOW   RATE   AND  ACCUM.    INFLOW   TO  POINT  2 
QI2=QI+QP2 
IF(OI2)170*187,187 
17C  QI2=.0 

187  DI= (QI2+QI 1 ) /2.*DELT 
ADV=ADV+DI 

DI=(DI/12.1>*12./ ( AREA-A2 ) 
A  I  I = A  I I+DI 

C   WRITE   POINT   2  OUTPUT 
012=01 2+. 005 
G2=G2+.001 
PI=PI+.001 
QP2=QP2+.005 

WRITE    ( 3  »92 ) NR  » NWS  *  NY  *  NMO ♦ NDAY  »T2»G2»PI » VAF2  * 0P2 » AC PS  »Q I  2  » A  I  I 

C   ADVANCE   POINT   1    TO  POINT   2   AND   ACCUMULATE   STORM  PRECIP. 
NER=NER+1 
G1=G2-.001 
RF1=BF2 
QPl=QP2-.005 
VAF1=VAF2 
T1  =  T2 
PA1=PA2 
01  1=012-. 005 

IF    (T2-24.0)  210.199*210 

199  JD=2 

C   COMPUTE   DAILY   SUMMARY    FOR    SUBDIVIDED   DAYS   AND  WRITE   ON  DISK. 

200  G=G2+.C01 
ADV=ADV/24. 

Z (NMO*NDAY) =ADV 

WRI TE ( 2 ' UOB ) NMO ♦ NDAY  *G ♦ ACP ♦ VAF2  » AQP T . AP V ♦ ADV  » NUM 

NUM=NUM+1 

GO   TO    (  110*206)  » JD 

206  IF( IMO ) 207*400*207 

C   REZERO  ACCUMULATORS 

207  ACP=0.0 
APV= .0 
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ADV=.0 
AQPT=.0 
T  1  =  •  0 

C  ADVANCE   POINT  2   TO  POINT  3 
210  VAF2=VAF3 
A2  =  A3 
G2  =  G3 
BF2=BF3 
T2  =  T3 
PA2=PA3 
P2  =  P3 
NMO= I  MO 
NDAY= I  DAY 

IF(T2-24. )133tl33»ll6 

C  SUMRE   PRINTS   DAILY  MONTHLY   AND   ANNUAL  SUMMARIES 
<*00  NMO=0 

WRI TE ( 2 1  JOB ) NMO »NDAY »G . ACP  » VAF2  »AOPT  * APV  »ADV  »NUM 
CALL  SUMRE( T ♦ Z » NR » NWS ♦ NY ♦ GHS ♦ VM I N ♦ Y A » YV ) 

GO  TO  9  5 
END 
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APPENDIX  F.-FORTRAN  LISTING  FOR  SUBROUTINE  CHECK 


SUBROUT INE   CHECK ( Y V » YQ ♦ NR  *  NWS ♦ NY  » VM I N * QM I N  » BP  *  BE  *  ARE A  » NUM  » Y A  * GHS ) 

C  READ  AND  STORE   INPUT   ON  DISK 

DIMENSION YV (40)  . YQ ( 1 00  )  ♦  YA ( 40 ) 
COMMON   JAB  *  JOS 

1  FORMAT ( 8F8.  2 ) 

2  FORMAT (I3.I4.I3.4F6.2.F8.1.I4) 

3  FORMAT ( 10X»2I3rF5.1»3F6.2) 

4  FORMAT ( 1H  8F8  .2 ) 

5  FORMAT (  1H1 > 

6  FORMAT ( 1HO    ' VOLUME   RATING   AT  ONE   FOOT    INTERVALS    IN  ACRE   FEET  • ) 

7  FORMAT ( 1 HO    'PRINCIPAL   SPILLWAY   RATING  AT   0.2   FOOT    INTERVALS   IN  CF 
IS  '  ) 

8  FORMAT ( 1H0' SURFACE   AREA  RATING  AT   ONE   FOOT    INTERVALS    IN  ACRES') 

10  WRITE (3.5) 
H  =  1 

K  =  8 
ZA  =  0 

WRITE(3»6) 

1 1  READ (2»1)  ( Y V (  I)  » I  =  N  »K ) 
WRITE ( 3  »4 !  ( YV ( I  )  » I=N  »K ) 
IF(YV(;<!  )  17. 18  1 17 

17  K=K+8 
N|  =  N+8 

GO  TO  11 

18  N  =  l 
K  =  8 

WRITE( 3.7) 

20  READ(2.1) ( YQ ( I) »I=N»<) 
WRITE(  3.4)  (  YQ  (  I  )  .  I  =N  .  K.  ) 
IF(YO(K)  )21»22»21 

21  K=K+8 
N  =  N  +  8 

GO  TO  20 

22  N  =  l 
K  =  8 

WRITE (3.  8) 

23  READ ( 2  . 1 )  ( YA(  I )  ♦  I=N»K ) 
WRITE( 3.4)  ( YA ( I ) » I =N  »K ) 
IF( YA(<) 124.30.24 

24  K=K+8 
N  =  N+8 

GO  TO  23 

30   READ ( 2  »  2 )NR .NWS.NY  » VM I N  »QM IN .BP .BE  »  AREA .NUM 
JAB=1 

3  5   R  EAD ( 2 » 3 ) NMO » NDAY » T I M ♦ GAGE  » P . BF 

WRITE ( 1 1  JAB  JNMO.NDA Y.TIM .GAGE. P.BF 
IF ( NMO )40. 40.35 
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40  JAB=1 

GC   TC  70 

CHECK    INPUT   DATA   FOR   TIME   SEQUENCE  AND  GROSS  GAGE. HEIGHT  ERROR 

51  FORMAT (1H   I  3 » 14 ♦ 13 » 2  I  3 ♦ F6 . 1 » 3F 6 . 2 ) 

52  FORMAT ( 1H05 X , 44HGAGE   HEIGHT    INCREMENT   GREATER   THAN   0.40  FEET) 

53  FORMAT ( 1H05X » 45HCARDS  MISSING  OR   DATE-TIME  GROUP  OUT   OF  ORDER) 

54  FORMAT (1H0) 

55  FORMAT ( lH   I  3  •  1 4 , I  3 » 2  I  3  ,  F6 .  1  * 3F6 • 2  ) 

56  FORMAT    (lHO'LAST   DAY   NUMBER  IS'»I4) 
65  ICO=ICO+l 

6  8  w R I T E ( 3 ♦ 5  5 ) NR ♦ NW S  »  N  Y ♦ L MO  *  L  DA Y  »  T I M  » G  *  P  »  BF 1 

WRITE(3»51)NR  »NWS  *  NY ♦ NMO  » NDAY  »  T I  ME ♦ G2  *  P 2  » BF2 

WRITEJ 3t54) 

GO  TO  72 
70  ICO=NUM-l 

72   READ ( 1 • JAB) LMO » L DA Y » T I M ♦ G 1 »P1 »BF1 

IF ( LMO ) 74 , 118*74 
74    IF    ( TIM-24.0 > 97,93 »80 

76  LMO*NMO 

L  DAY  =  1 

GO  TO  80 
78  LDAY=NDAY 

T  I M  =  T  I  M  E 
80  ICO=ICO+l 

32   READ ( 1 • JAB > NMO, NDAY* TIME, G2»P2 »BF2 

IF(NM0)84,118»84 
84  Dir=Gl-G2 

P  =  P1 

G  =  G1 

G1  =  G2 

P1  =  P2 

IF(ABSIDIF)-. 40)86,86,85 

8  5    (  I F ( T I  ME- ( T I M+ • 1 )  ) 9 1 » 8  6  »  9 1 

86  IF(TIME-24.0) 108,93.87 

87  IF(NDAY-LDAY-1)89,88»93 

88  IF(NM0-LM0)93 #78,93 

89  IF(NMO-LMC-1)93,90»93 

90  IFtNDAY-l )93,76,93 

91  WRITF(3»52) 
ZA=1 

IF( TIME)92»65  »92 

92  IF1TIME-25. 0)68,65,68 

93  WRITE(3»53) 
ZA=1 

GO  TO  68 

94  LMO=NMO 
LDAY= 1 
ICO=ICO+l 
GO  TO  96 

95  LDAY=NDAY 
ICO= ICO+I 

96  T IM=T I  ME 

9  7  R FAD ( 1 • JAB) NMO, NDAY, TIME, G 2, P2  »BF2 

IF( NMO)  9 3,1 18, 90 
98  DIF=G1-G2 
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G  =  G1 

P-Pl 

G1-G2 

P1=P2 

RF1=BF2 

IF(ABS(DIF)-.40)100*100*99 
99    IF(TIME-(TIM+«1) )91»100*9l 
100    IF( TIME-24.0) 102*93.114 
102    IF ( T IME-TIM ) 1 10*93 »  104 
104    IF(NDAY-LDAY) 93*106*93 
106    I F ( NMO-LMO > 93 »96 »93 
108    ICO= ICO-1 

IF(TIME)93»93tllO 
110  IF(NDAY-LDAY-1)112»116»93 

112  IF ( NMO-LMO- 1 > 93*113*93 

113  IF(TIME)93»94»93 

114  IF(i\DAY-LDAY)93*115»93 
lib    IF( TIM193, 68*93 

116  IF(NMO-LMO)93»95»93 

118  WRITE( 3*56)  ICO 
IF(ZA) 119*120*119 

119  CALL  EXIT 

120  RETURN 
END 
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APPENDIX  G.-FORTRAN  LISTING  FOR  SUBROUTINE  SPPNT 


SUBROUTINE   SPPNT ( Y tN »X  ♦  Y  I  ) 

SIX   POINT   PARABOLIC    INTERPOLATION     COEFFICIENTS   ARE   NOT  STORED 
VERSION    II    DISK   MONITOR    IBM  1130 

Y    IS   AN  ARRAY  WHICH   CONTAINS   DATA   POINTS  AT   EQUALLY   SPACED  INTERVALS. 
BETWEEN  WHICH    INTERPOLATION    IS   DESIRED.      Y ( N )    IS   THE   VALUE   WHICH  IS 
AT    THE   LEFT   OF   THE    INTERVAL    IN   WHICH    INTERPOLATION   WILL   BE  MADE. 
X    IS   THE   FRACTIONAL   DISTANCE   TO  WHERE    INTERPOLATION    IS  DESIRED* 
It    E.      IF    INTERPOLATION    IS   DESIRED  AT   A   POINT   HALFWAY  BETWEEN 
Y(N)    AND   Y(N+1).   X   WOULD   EQUAL    .5.      Yl    IS   THE    INTERPOLATED  VALUE 
WHICH    IS  RETURNED. 
DIMENSION  C(6»6)  »A!6)  »Y(  1  ) 

DATA   C/432.  »-9l8.  »765.  »-313.  »63*  »-2'040«  »4436«  »-37  54»  »1551.  »-31 

14. *2  5. »4C80»  *-8752 . »7414 • »-3078. *626»  »-50« »-4080. » 871 2 • »-73  56 • t 305 
28 • *-624. »50«  »2040«  »-4346. »366l • »-l521. »311. »-25»  »-408 • f 868*  »-730. ♦ 
3303. .-62. »5./ 
IriN-2>  50.50.70 
50  WRITE  (3.60) 

6D   FORMAT ( 1 0   ERROR    IN  SPPNT'! 

CALL  EXIT 
70  YI=0. 

00  9  0   J =1.6 

A  (  J  !  =  C  (  6  .  J  ! 

DO  80   <  =  2.6 

A( J) =A ( J)*( 2.+X) 

L  =  7-K 
80  A ( J ) =A ( J )+C ( L  .  J ) 

M=N-3+J 
9G   Y I =Y I +Y ( M ) *A ( J ) 

YI=YI/24. 

RETURN 

END 
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APPENDIX  H.-FORTRAN  LISTING  FOR  SUBROUTINE  SUMRE 


PR  SP  RAIN 
PR   SP  RAIN 
(ACFT)  (ACFT) 


SUBROUT  INE   SUMRE  (  T  ♦  Z  »  NR  *  NWS  t.NY  »GHS  »  VM I N  »  YA  t  Y V ) 
DIMENSION  Z ( 12»31 )  »T ( 12  ) 
COMMON   JAB » JOB 

1  FORMAT ( 1H 1  ) 

2  FORMAT ( 1H050X »' TOTAL    I NF LOW= 1  . F 1 0 • 2 t IX  *  1 SFD 1  ) 

3  FORMAT ( 1H    1    NR  STA  YR  MO  DY     AREA       GH     RAIN  VOL 
1  I  NFL   N'JM  LOSS  1  > 

4  FORMAT ( 1H 1  *    NR   STA   YR   MO  DY     AREA        GH     RAIN  VOL 
1  INFL   NUM  LOSS  1  ) 

5  FORMAT ( 1H    '  (ACRES )    (FT)    (IN)  (AC^T 
1                                  (ACFT  )  (ACFT  )    (ACFT)        (  IN)  1  ) 

6  FORMAT (1H   I  3  »  I 4 » 3  I  3  » F6 • 2  » F6 • 2 » F 5 • 2 » F8 • 2  » F6 • 2  » F7 • 2  » 14X  » F7 • 2  » 1 4 » 5X » 
13F7.2 ) 

7  FORMAT    (1H  »12F10.2) 

8  FORMAT ( 1H  F 1 0 . 2 ♦ 1 0 X » 1 0 F 1 0 . 2  ) 

9  FORMAT    (1h   » F 1 0 . 2 . 1  OX . F 1 0 . 2 ♦ 1 0 X » F 1 0 . 2 . 1  OX . F 1 0 . 2 » F 1 0 . 2  *  1  OX » F 1 0 . 2 » 1 0 
1X.F10.2  ) 

10  FORMAT  (1HJ) 

11  FORMAT    (1H0120M  JAN  FEB  MAR  APR  MAY 
1            JUNE              JULY                AUG              SEPT                OCT  NOV 
2DEC) 

12  FORMAT  (  1H  4  OX  ♦  1  S  I  T  E  '  *  I  4  ,  •      WATERSHED  1  *  I  5  »  1    INFLOW    (SFD)    IN  19SI2) 

13  FORMAT  (1H    1  " 

!   ) 

14  FORMAT (1H  16X.F6.2»6X»F5.2»F8.2»F6.2.F7.2»14X.F7.2.I4»5X.F7.2»7X» 
1F7.2  ) 

15  FORMAT (  1H0 ' ANNUAL  SUMMARY • 1 14X »F5 .2 ♦ F3 • 2 »  F6 . 2 » F7  .  2  ♦  14X  ,  F 7 . 2  ♦ 
19X,F7.2»7X»F7.2) 


C  INITIALIZE  ACCUMULATORS 
AA1=.0 
AA2=  .0 
AA3=.0 
AA4= .0 
AA5=.0 
AA6=  .0 
AA7=  .0 
NMO  =  0 

NAD  =  GHS-VM.IN 
XAD=GHS-FLOAT ( IFIX (GHS ) ) 
CALL   SPPNT ( YV »NAD*XAD » VI ) 
ACUM  =  .  0 
N  =  0 

ARAIN=.0 
APRSP= .0 
ARN=.0 
AINFL=.0 
C=24/12. 1 
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A  A  F  L  =  •  0 
AETSI=.0 
A  F  L  M  =  •  0 
A AFL  =  •  0 
ETSI=.0 
W  R I  T  E  (  3  »  4 ) 
WRITE (3. 5) 
J0B=1 

READ ( 2 • JOB) NM»ND*GA.ACP  »VAF2  * AGPT »APV»ADV»NUM 

M=IMM 

V1=VAF2 

VOLl=VAF2 

GO   TO  30 

25  READ(2*J0B)    NM ♦ ND » GA » ACP ♦ VAF2 * AOPT ♦ APV ♦ ADV *NUM 

IF(NM-M)40»30»40 
2  9  WRITE  C  3»3) 

to  R I  T  E  (  3  »  5  ) 

C   COMPUTE   AND   PRINT   DAILY  SUMMARY 
30   GHA= (GHS+GA ) /2. 
NAD=GHA-VMIN 
XAD  =  GHA-FLOAT ( IF  IX (GHA  )  ) 
CALL  SPpNT(YA»NAD»XAD»AR) 
GHS=GA 
ADV= (ADV*C) 

AFL=APV+ADV-AQPT-( VAF2-V1 ) 

ETS   = ( AFL*12. /AR ) 

ETSI=ETSI+ETS 

AAFL=AAFL+AFL 

AFLM=AFLM+AFL 

V1=VAF2 

ACUM=ACUM+AR 

N  =  N+1 

ARAIN=ARAIN+ACP 

VOL2=VAF2 

APRSP=APRSP+AQPT 

ARN=ARN+APV 

AINFL=AINFL+ADV 

AETSI=AETSI+ETS 

ACP=ACP+.001 

VAF2=VAF2+.005 

AQPT=AQPT+.005 

APV=APV+.005 

ADV=ADV+.005 

WRITE! 3  »6 )NR  *NWS  »NY»NM  »ND  » AR  »GA  »ACP  » VAF2  »AQPT  » APv ♦ ADV  »NUM  * 
lAFL.AAFLtETSI 

GO  TO  25 

C  COMPUTE   AND   PRINT   MONTHLY   SUMMARY  COMPUTATION 
40  WPI TE ( 3  » 13 ) 

DVOL=VOL2-VOLl 
AN=N 

ACUM=ACUM/AN 
AA1=AA1+ARA IN 
AA2=AA2+DVOL 
AA3=AA3+APRSP 
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AA4= AA4+ARN 
AA5=AA5+AINFL 
AA6= AA6+AFLM 
AA7=AA7+AETSI 
ARAIN=ARAIN+.001 

WRITE ( 3  ♦  14)ACUM»ARAIN»DV0L*APRSP»ARN»AINFL » N ♦ AF LM » AET S I 

WRITE(3»13) 

VGL1=V0L2 

M  =  M  M 

ACUM=.0 

N-0 

ARAIN=.0 

APRSP=.0 

ARN=.0 

AINFL=.0 

AFLM= • 0 

AETSI=.0 

IF(NM) 50*50*2  9 

C   COMPUTE   AND   PRINT   ANNUAL  SUMMARY 

5  0  WRITE! 3 » 15 ) AA1 »AA2  *AA3  » AA4 » AA5  ♦  AA6 ♦ AA7 
WRITE ( 3.1 ) 

WRITE(3»12)NR»NWS*NY 

WRITE ( 3*11 ) 

DO  52  NF=1,12 

DO   5  2  M=l,31 
52  Z(NF»M)=Z(NF»M)+.005 

M-l 

NF  =  1 
54  K  =  0 

56  WRITE    (3*7) (Z(NF*M) »NF=1»12) 
M  =  M+1 
K  =  K.+  1 

IF  (M-29)58*62»58 
58   IFCK-5 ) 56*60*56 
60  WRITE  (3*10) 

GO  TO  54 

62  WRITE    ( 3  •  8 ) Z ( 1*29)  » ( Z ( NF  » M )  *NF  =  3»12) 
M  =  M+1 

WRITE  (3»8)Z(1»30)»(Z(NF.M)»NF=3*12) 
M  =  M+1 

WRITE    ( 3»9)Z( 1 »31 ) *Z ( 3  *31 ) *Z( 5  »31 ) »Z (7*31 ) »Z { 8  »31 ) »Z 1X0*31 ) »Z( 12 
1,31) 

SUM1=#0 

TSUM=.0 

DO  66  NF=1»12 

DO  64  M=l*31 

TSUM=TSUM+Z (NF»M)-.005 
64  SUM1=SUM1+Z (NF»M)-.005 

T  (  NF  )  =SU,Ml+«005 
66  SUM1=«0 

WRITE    (3*7) (T(NF) »NF= 1 » 12 ) 

WRITE( 3»2)TSUM 

RETURN 

END 
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